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A GENERALIZED POWER LINDLEY DISTRIBUTION WITH 
APPLICATIONS 


GAYAN WARAHENA-LIYANAGE AND MAVIS PARARAI 


ABSTRACT. The Exponentiated Power Lindley (EPL) distribution which is an extension 
of the Power Lindley Distribution is introduced and its properties are explored. This new 
distribution represents a more flexible model for the lifetime data. Some statistical properties 
of the proposed distribution including the shapes of the density and hazard rate functions, 
the moments and moment generating function, skewness and kurtosis are explored. Entropy 
measures and the distribution of the order statistics are given. The maximum likelihood 
estimation technique is used to estimate the model parameters and finally an application 
of the model with a real data set is presented for the illustration of the usefulness of the 
proposed distribution. 


1. INTRODUCTION 


The power Lindley (PL) distribution was proposed by Ghitany et al. [4]. This distribu- 
tion is an extension of the Lindley (L) distribution which was proposed by Lindley [7] in 
the context of fiducial and Bayesian statistics. Properties and applications of the Lindley 
distribution have been studied in the context of reliability analysis by Ghitany et al. [3]. 
Several other authors including Sankaran [13], Asgharzadeh et al. [1] and Nadarajah et al. [8] 
proposed and developed the mathematical properties of the generalized Lindley distribution. 
The probability density function (pdf) of the Lindley distribution is given by, 


(1.1) uD = 


Using the transformation X = Ya, Ghittany et al. [4] derived the power Lindley (PL) 
distribution given by 


(lt+yje, y>0,8>0. 


D 


2 
(1.2) f(x;a, = = z0 Ee gle P, gS 0, > 0/8 > 0. 

The survival function and cumulative distribution function (cdf) of the power Lindley dis- 
tribution are 


(1.3) S(x) = (1 + =) a 
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and 


(1.4) Fac ea eee Co ele (: i. FS em 





for x > 0,a, 8 > 0, respectively. 

The purpose of this paper is to develop a three-parameter alternative to several lifetime 
distributions including the gamma, Weibull, exponentiated Weibull, exponentiated Lindley, 
and lognormal distributions. In this context, we propose and develop the statistical proper- 
ties of the exponentiated power Lindley (EPL) distribution and show that it is a far better 
model for reliability analysis. 

Our aim in this paper is to discuss some important statistical properties of the EPL 
distribution. This discussion includes the shapes of the density function and hazard rate 
function, reversed hazard rate function, moments, moment generating function, distribution 
of order statistics and model parameter estimation by using the maximum likelihood method. 
Finally, applications of the model to real data sets in order to illustrate the applicability and 
usefulness of the EPL distribution are presented. 

This paper is organized as follows: In section 2, the model and some of its statistical 
properties including shapes and behavior of the hazard function are presented. Distribution 
of order statistics, Moments and related measures are given in section 3. Section 4 contains 
entropy measures. In section 5, mean deviations, Lorenz and Bonferroni curves are presented. 
In section 6, we present the maximum likelihood method for estimating the parameters of 
the distribution. Applications are given in section 7 and thereafter the concluding remarks. 


2. THE MODEL, SUB-MODELS AND PROPERTIES 


In this section, the cdf, pdf, hazard and reverse hazard functions, and their shapes are 
presented. We will define G(x) = |[F(ax)|* for w > 0 where F(x) is the cdf of the power 
Lindley distribution given in (1.4). 

(2.1) G(2;a,8,w) = h = (14 = - 


We refer to G(x) as cdf of the EPL distribution. Plots of the cdf for the EPL distribution 
for several values of the parameters, a, 8 and w are given in Figure 2.1. 





j , z>0,a>0,8>0,w >Q. 
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FIGURE 2.1. Plot of the CDF for different values of a, 8 and w 


The pdf of the EPL distribution is given by 


(z;a,ßb,w) = CBW el L gegote Po 1—{1+ ea en E 
g a Bale = B L 1 B a 1 , 
for x > 0,a > 0,8 > 0,w > 0. Plots for the pdf of EPL distribution are given below and 
each plot has been generated by fixing one parameter at a time. 
e a is fixed: see figure 2.1. 
e £ is fixed: see figure 2.2. 
e w is fixed: see figure 2.3. 
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FIGURE 2.3. Plot of the PDF for different values of a,w and 8 = 2.0 
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FIGURE 2.4. Plot of the PDF for different values of a, 8 and w = 1.8 


The pdf of the EPL distribution is unimodal. It increases and decreases for various values 
of the parameters giving the shapes obtained in the above plots. The plots in figs 2.3 and 
2.4 show that for a < 1, the pdf is decreasing and for values of a > 1, the pdf is unimodal. 
For values of w > 1,8 > 1, the graph seems almost symmmetric. 


2.1. Sub-models. Some sub-models of the EPL distribution are presented in this section. 
e When w = 1, we obtain the PL distribution. 
e When a = 1, we obtain the exponentiated Lindley (EL) distribution. 
e When w = a = 1, we obtain the Lindley (L) distribution. 
e When (6 = 1, we obtain exponentiated two-component mixture of Weibull distribu- 


tion (with shape parameter a and scale 1) and a gamma distribution (with shape 
parameter 2a and scale 1). 


e When a = 2, we obtain an exponentiated two-component mixture of Rayleigh distri- 
bution (with scale 8 and scale 2) and a gamma distribution (with shape 4 and scale 


8). 


2.2. Shapes and Stochastic Orders. In this section, we present the mode and discuss the 
shape, as well as stochastic orders of the EPL distribution. To obtain the mode, we solve 
the equation ace) = 0, for x. Note that 


log(g(z)) = log(a) + 2log(Z) + log(w) — log(1 + 8) + log(1 + x°“) — Bx 
(2.2) + (a—1)log(x) + (w —1)log h — (1 + Ee ee 





B+1 
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aint antala) 


so that = 0, results in 


a-1 i 2 a) pa—-1,—Bx™ 
aT abrt? Leisi ap (1 + 2%)x2%*e = 
1+a° a (1+ 6) — [1+ 8(1 + 2?) ef 

Note that lim, g(x) = oo, and lim, 4. g(x) = 0. 

Let X; be distributed according to EPL(a, 8,w), with cdf and pdf G; and g;, respectively 
for i = 1,2. We say X; is stochastically greater than X; in likelihood ratio if go(x)/gi(x) is 
an increasing function of x. It is well known that likelihood ratio order implies failure rate 
order which in turn implies stochastic order, see Shaked and Shanthikumar [12] for additional 
details. 

e If 6, = Py and a, = ag, then Xə is stochastically greater than X, with respect to 
likelihood ratio order if and only if w2 > wy. 
e If a, = ae and wı = w2 then Xə is stochastically larger than X, with respect to 
likelihood ratio order if and only if 6, > do. 
Note that 


= 0. 





g2(x) (1+ B1)arg?wofo(1 + 2°?) 702701 er part 
gi(x) (1 + b2)a1?w b11 + £1) 


h = (1 aie Box? eon] 
(2.3) x Ted- 


Bir™ at wi—1' 
Ba aie 


If 61 = Bo, and ay = ag, then 


(2.4) Kae = h 7 (1 4. ET) ep] 























and is such that 
wlw = Ww) 


K(x) = - | sa oa" exp(—B0*) 


Bx Aa w2—w1—1 
‘earn ia Be >0 
| ( “Beis aac 


2.3. Quantile Function. The quantile function is the solution of the equation 
Í z (- p+ px 

B+1 
Thus, the quantile function, say Q(p), defined by G(Q(p)) = p is the root of the equation, 


h 2 (- PLF eer) exp(-80(0)")] EA Oecd: 





x 





(2.5) 


if and only if wo — w1 > 0. 


Jer] =p where0<p< 1. 





6+1 
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Let Z(p) = —1 — 2 — BQ(p)*. 


We have 
j+ (52) etzo ESES] ita i 
1+ (74) exp{Z(p) 41+ 6} = pil 
Z(p)exp{Z(p)} = ~(6+10-p 





exp(1 + 8) 
So the solution for Z(p) is 

Z(p) =W (—(8 + 1)(1 — p°) exp(-1 — 6). 
for 0 < p < 1, where W(.) is the Lambert W function, [2]. 


Now, 
-1—B-BQ(p)* = W (-(8+1)(1-— p“) exp(-1 — 8)) 
BQ(p)* = —1-8-W (-(6 +1)(1-— p”) exp(—1 — §)). 
So we have 
(2.6) Q(p) = e m E ne me) c-e) 
: p) = 3B p“) exp : 


2.4. Hazard and Reverse Hazard Functions. The survival function for the EPL distri- 
bution is given by, 


G(z;a,8,w) = 1-G(z;a, B,w) 


px“ — pre z 
i-[1- (14+ Je | ; 


The hazard and reverse hazard functions are given by 





(2.7) 











àa(x; a, b, w) = Gee) 
G(x; a, B,w) 
2w a) »a—1,— px —Ba ee 
re ea ten h- (14+ 8 )e f T 
= Bx® oe wW ’ 
1 h 0+ e | 
and 
l _ g(zr;a,ß,w) 
a(x; a, p, ) m OG a, B,w ) 
a “(14+ 2%)xete be" £ = (1 + gz) Em 7 











Bae E w , 
= (e | 


respectively. Plots of the hazard function are given below: 
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FIGURE 2.5. Plot of the hazard function for different values of a, 8 and w 


The graphs of the hazard function for four combinations of the values of the model pa- 
rameters show various shapes including monotonically increasing, monotonically decreasing, 
uni-modal, bathtub, and upside down bathtub shapes with four combinations of the values 
of the parameters. This attractive flexibility makes the EPL hazard rate function useful and 
suitable for non-monotone empirical hazard behaviors which are more likely to be encoun- 
tered or observed in real life situations. 


3. MOMENTS, MOMENT GENERATING FUNCTION AND RELATED MEASURES 


In this section, moments and related measures including coefficients of variation, skewness 
and kurtosis are presented. A table of values for mean,standard deviation, coefficient of 
variation (CV), coefficient of skewness (CS) and coefficient of kurtosis (CK) is also presented. 


3.1. Moments. The rt moment about the origin of a continuous random variable X, de- 
noted by w., is , 


w. = E(X”) = f z”glxjdx forr =0,1,2,... 
In order to find the moments, consider the following lemma. 
Lemma 1 

Let, 





E a w—1 
Ly (a, B,w,r) = / x”(1 + a ls h — (1 + = z) oe e-l% dr 


ee) Qa w—l1 
f (1 + ojo f n (1 ft = z) = e= dr, 
0 
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then, 
SE fw -1) fi) (G41) (-D BT Rt rat t+D 
Li(a, 8,w,r) = oe) | = EN 
(a, 8, w,r) dada ds i y k J aB EIGE DE 


(Pl a-y ( ; ‘) is, 
we have 


Lila, b,w,r) = D (* 7 J (—1)' HE] ee fo + reret le p dy 





= D (07 ER Po teen pa Ha ean 
i 0 

= = w—1 (l : i fS a\j+l,atr—1 ,(-iBr*—Bx%) 

= > i aA, R E a 
Tn w-—1 a j+1 (—1) 8 E atak+r—1,—(i+1) Ba 

ge a ra ae 


Now consider, 


(3.2) i. getak+r=1 o= (G41) Bae de. 
0 











1/a 
tastei then @ a e aa | |. 
et u= p(i+ 1)x en & =ap(i + 1)x™t and x n : 


The above integral can rewritten by using the complete gamma function T (a) = f xole tda 











as, 

a (RE) au o de e u 

0 L(BG+1) apli + 1) natan 
= T(k+ra™!t +1) 
T a 
Consequently, 

LAE eN (8) tN ETE ra! +1) 
$ = . 
i(a, B,w, n) 2 ar E ( 4 ) (C) ( k o + 1) [B(i + POREN 














Now using Lemma 1, the r? moment of the EPL distribution is given by 
aB?w 


(3.3) Hi, = EX) = S 


Lila, B,w,7). 
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The first four moments of X can be written as, 

















2 

m = BCX) = Ello gw) 
2 

Hy = BOX?) = FAI (a, 8,0,2). 
2, 

m = BUX) =F Ad (a, 8,03). 
2 

Hy = BUX!) = SEE Tilo 8,04). 


The mean, variance, coefficient of variation (CV), coefficient of skewness (CS) and coefficient 
of kurtosis (CK) are given by 





2 
(8.4) p= m = BCX) = FIs a, pw, 1), 
(3.5) o?’ = Ha _ pr, 





TE ,,2 1 
(3.6) EE E 
u u u 


(3.7) cg- El = )") _ os = 3um + 2H 


[E(X — p)” (Hp — p?)’.? 





and 
(3.8) cg- E- n) = Aups + 6p? u — one 

[E(X — p)?] (u3 — p?) 
respectively. Table (3.1) presents the mean, standard deviation (SD), coefficient of variation 
(CV), coefficient of skewness (CS) and coefficient of kurtosis (CK) for some values of the 
parameters a, 8 and w. Note that CV, CS and CK do not depend on the parameter a. 





GENERALIZED POWER LINDLEY DISTRIBUTION 





Mean 


SD 


CV 


CS 


CK 





2.0000 
9.8765 
80.6584 
0.2880 
0.4800 
1.3548 
0.0875 
0.0714 
0.1000 


0.8165 
NaN 
NaN 

0.4910 

0.7494 

1.0488 

0.2377 

0.2688 

0.3644 


1.1547 
1.5411 
1.8439 
0.3536 
0.8062 
1.0801 
NaN 
0.2887 
0.6547 


1.5651 
NaN 
NaN 

1.3567 

1.8670 

0.3672 

1.7858 

2.0456 

1.9691 


2.2913 
NaN 
NaN 

2.5053 

2.0369 

1.5013 

3.7779 

4.6601 

4.2695 





1.4361 
5.7117 
40.2338 
0.2991 
0.3978 
0.9641 
0.1127 
0.0731 
0.0877 


NaN 

NaN 

NaN 
0.4456 
0.5672 
0.6522 
0.2734 
0.2570 
0.3038 


0.4668 
0.6478 
0.7569 
NaN 
0.3409 
0.4907 
NaN 
NaN 
0.2816 


NaN 
NaN 
NaN 
0.9267 
0.8657 
NaN 
1.4336 
1.8063 
1.7887 


NaN 

NaN 

NaN 
1.3973 
1.3021 
1.3411 
2.3570 
3.4579 
3.3849 





1.2192 
4.2016 
26.9248 
0.3174 
0.3634 
0.7988 
0.1363 
0.0759 
0.0824 











NCIWN OWN OWIN OWN OWN OWN OWN OWN OWE 





NaN 

NaN 

NaN 
0.4509 
0.5032 
0.5159 
0.3180 
0.2610 
0.2810 





0.3345 
0.4742 
0.5525 
NaN 
0.2450 
0.3633 
NaN 
NaN 
0.2026 





NaN 
NaN 
NaN 
0.8706 
0.8288 
NaN 
1.3956 
1.7837 
1.7724 





NaN 

NaN 

NaN 
1.2877 
1.2339 
1.3525 
2.2156 
3.3524 
3.3090 











NaN: Not Defined 


TABLE 3.1. Table of Mean, SD, Coefficient of Variation, Skewness and Kurtosis 
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3.2. Conditional Moments. For lifetime models, it may be useful to know about the 
conditional moments which can be defined as E(X” | X > x). In order to calculate these, 


we consider the following lemma: 
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Lemma 2 
Let, 
a r ay) ,,a-1 px“ — Br” a — Bx 
Lala, b,w,r,t) = z (1+ °)z 1—(1+ e e dx 
t +1 


lo) Q w—l1 
1 a),atr-1 es 1 Bx —Ba~ —Ba~ : 
/ (l+2%)z | ( Pai e e” dg 
then, 


rosa IA OE e a 


i=0 j=0 k=0 











Proof. Using the same procedure that was used in Lemma 1, this can be simplified into the 
following form. 


i j+l . S 
Tilo esos > ee — ') ( ) ( i ') 5 2 + L gatak+r-1e-(i+1)b1® Jy. 


Now consider, 





oo 
— aig Q 
f gataktr le (i+1)8x dx. 
t 


1/a 
Let u = (i+ 1)x%, then & = aB(i+1)x° and z = z] l 


The above integral can be rewritten by using the complementary incomplete gamma function 
T (a,t) =f. redz as, 


m | u | (a) ou du E i yktrea emu i, 
piya LAG + 1) aB(i + 1) ater a [B(t + Lee 
Tik+ra7*+1,8@+ 1)t) 


a a 











Consequently, 


aesan SEE MCT NEE Ee 


i=0 j=0 k=0 

















Now using Lemma 2, the r“’conditional moment of the EPL distribution is given by 
ap?’w L(a, p, wW, T, x) 








PAE grea A 
aa. oP wha, Bwr a) - 1 
7 1 Bar —Bar& 
p+ (1+ gee B 


a3?w Lala, B,w,r, £) 
B+10 +8 pr) AA 
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The first four conditional moments are given by, 
aB?w Lo(a, B,w, 1,2) 


























WANK >a) = oa a pro) e Ar 
BX >a) = Se pee 
BX >n = SS aeit 
BAX >a) = e he 














The mean residual lifetime function is given by E(X|X > x) — z. 


3.3. Moment Generating Function. The moment generating function (MGF) of a con- 
tinuous random variable X, where it exists, is given by, 


OO 


(3.9) Mx(t) = E(e*) = / e'g(x)dz. 


OO 


The MGF of the EPL distribution is Da by, 


(3.10) = - oat “ L(a, bwn). 





Proof. Consider, 


2 oo a w—l 
Mx(t) = a | eta) ete h — (: + = z) = dx. 


This can be simplified into, 








Bw oo i J+ 
NM) = geet oe 


i=0 j=0 k 


COC ae 


oo 
x f gotok-le— (+1) be et dx. 
0 








(3.11) 


The proof of equation (3.10) is similar to the proof of Lemma 1, but without using the 
definition of the gamma function. Now consider, 


oo 

— —(4 a 
| getak le (i+1)8x edz. 
0 





OO 
: : : le 
Using the series expansion, e’” = ys , we have 
n=0 
f getek- 1 pC iBa%—Ba) eae = [S ye x” getek- 1 e iBx% —Ba® dz. 
nl 
0 n! 


Consequently, the MGF of the EPL distribution ae to 


D aß? oo tr 
(3.12) Mx(t) = reel > Lala, p, w, n). 
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3.4. Distribution of Order Statistics. Order Statistics play a vital role in probability 
and statistics. In this section, we present the distribution of the order statistics for the EPL 
distribution. The pdf of the i order statistic is given by: 


gle) = oS) O - ata). 





Using the series expansion, 


we have: 








aß’ wn! ca an gre To (2-4 
SD De e D ') 


Lt Bt Bae’) pie hei ic 
xo h- (= Je | 


oe) 


- ee See (CE) 








rena Al O 
7 niea El 


D (SCE IOC 


42 pam+a-1 


(8 + Lert 





x (=t en Beh et) 


Now using the series expansion, 
2 BPP(k +1)? 


e Be (k+1) 
p! 


b] 


p=0 
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we have: 








9%) Git aa 2 5 2 > ("7") A 








j=0 k=0 l=0 m 
(—1)i tk gltp+2 gall+m p) 1( 
PERIE 





x 








4. MEAN DEVIATIONS, LORENZ AND BONFERRONI CURVES 


15 


In this section, we present the mean deviation about the mean, the mean deviation about 
the median, Lorenz and Bonferroni curves. Bonferroni and Lorenz curves are income in- 
equality measures that are also useful and applicable to other areas including reliability, 
demography, medicine and insurance. The mean deviation about the mean and mean devi- 


ation about the median are defined by 


Di) = [ Cetin 
and 


D(M) -[ | 2 — M | g(x)dz. 


respectively, where u = E(X) and M = Median(X) = G~'(1/2) is the median of G. These 


measures D(u) and D(M) can be calculated using the relationships: 





oo H 
(4.1) Di) = Glu) -2+2 | zglede = 246l) -2 | ag(w)ae 
H 0 
and 
oo M 
(4.2) D(M) = —u + 2 | xrg(x)dx = u — 2 | xg(x)dz. 
M 0 
Now using Lemma 2, we have 
2a87w 
Di) = 2uG(y) — 20+ EE halo, Bw, 1p) 
Bl 
and 
E 2a 8?w 
D(M) = per B+1 Lala, B,w,1, M). 
Lorenz and Bonferroni curves are given by 
Jo tg(t)at L(G(a)) 
4. L EA B = 
(4.3) (Ga) =F and BE) = = 
or 


q 1 q 
(4.4) Lip) = 7 f tg(t)dt, and Bl) = — I tg(t)dt, 
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respectively, where q = G~!(p). Now using the Lemma 2, we can re-write Lorenz and Bon- 
ferroni curves as 


1 q 
Blip) = — | tg(t)dt 
(p) nih (t) 


a + m mie a zglede 


2 
= Ee Ç — a Lol; 


and 


u 


: |S zoet- f zga) 


1 2 
= 7 Ç = SEE balo, Buu 1.0) a 


5. SOME MEASURES OF UNCERTAINTY 


In this section, we present Shannon entropy [{10],[11], as well as the Rényi entropy, [9] for 
the EPL distribution. The concept of entropy plays a vital role in information theory. The 
entropy of a random variable is defined in terms of its probability distribution and can be 
shown to be a good measure of randomness or uncertainty. 


5.1. Shannon Entropy. Shannon entropy is defined to be 


H [g(X;a, B,w)] = E|- log(g(X; a, ß6,w))] $ 
Thus we have 


H[g(Xsa,8,4)] = tog ER 

















B+ B? 
= og | a [Brite 8.0.0) 
T EN 
q=1 q 
iia So 
p 1 a, B,w, a) 
p=1 a=0 
+ unh 2 r 


(5.1) x Lala, Bwal +8); 
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5.2. Rényi Entropy. Rényi entropy is an extension of Shannon entropy. Rényi entropy is 
defined to be 


1 CO 
(5.2) Irg(v) = i log (/ [g(z; a, Bw)" , vÆl,v>0Q. 
Rényi entropy tends to Shannon entropy as v — 1. Note that by using the series expansion 
in equation (3.1), we have 





[vee = (ESD O O 


oo 
x | gratvatka ve Blj+v)x dr. 
0 











Now, let y = (j + v)x%, then 
oo i a a T(irtv+k— ut 
| grotvatka ve B(i+v)z dr = ( Ca) 
0 











Consequently, Rényi entropy is given by 





T(r +u+k-—(%)) | 
alag + oye Ce]! 

foro lw: 
6. MAXIMUM LIKELIHOOD ESTIMATION 


In this section, the maximum likelihood estimates of the EPL parameters a, 3 and w are 
presented. The log-likelihood of a single observation x of X from the EPL distribution is 


log(g(z)) = log(a) + 2log(B) + log(w) — log(1 + 8) + log(1 + z“) + (a — 1) log(z) 


























(6.1) — px” + (w—1)log h — (1 + = ; Jeepa). 
The partial derivatives of log(g(xz)) with respect to the parameters a, 8 and w are: 
log(g(x 1 x log(x 
BUCI L + loga) EG) ga log(a) 
L qo EEO (1 BEP Ga oele) 
bodes jen 
log(g(z)) _ 2 1 a 
OB 6 8+1 
N a= i) = 0+ Be) 


1-(1+ aoe 
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and 
o log(g(x)) 1 Ba? \ pra 
= + ] 1-{1 ; 
Ow wW 08 i 6+1 
The total log-likelihood based an a random sample £1, £2, ..-., Un, of size n is l = X`; log(g(x;)) = 


XL] 4 The elements of the score vector are: 


ol n € x? log(x;) 
ae ee) | sou 429 - os log (;) 


i=1 

















n — Bre Bx} log(xi oo _ Bae a 
(6.2) — (w—1) > ems : si — (1 Arie i (Bx? log(x;))) 
i=l 1— (1+ fie 





ol 2n n s 
a6 ~ p Bei” 











w Ti i pri L Baty opr? na 
(6.3) = ta = alee eg a aaa 
i=1 1- (1+ te — Bat 
and 
ol Bx? g 
6.4 ae E l 1—{1 — Bx; . 
ne oe s bbe GE | 


i=1 


The maximum likelihood estimates, Ô of © = (a, 3,w) are obtained by solving the nonlinear 
equations oe = 0; 2 = 0, and £ = 0. These equations are not in closed form and must be 
solved via iterative methods such as Newton-Raphson method. 
6.1. Asymptotic Confidence Intervals. In this section, we present the asymptotic con- 
fidence intervals for the parameters of the EPL distribution. The expectations in the Fisher 
Information Matrix (FIM) can be obtained numerically. Let Ô= (a, Ê, &)T be the maximum 
likelihood estimate of © = (a, 8,w)”. Under the usual regularity conditions and that the 
parameters are in the interior of the parameter space, but not on the boundary, we have: 
vn(ô®- 0) aN N3(0, I71(©)), where I(©) is the expected Fisher information matrix. The 
asymptotic behavior is still valid if [(@) is replaced by the observed information matrix 
evaluated at ©, that is J(@). The multivariate normal distribution N3(0,J(Q)~!), where 
the mean vector 0 = (0,0,0)", can be used to construct confidence intervals and confidence 
regions for the individual model parameters and for the survival and hazard rate functions. 
The likelihood ratio (LR) test can be used to compare the fit of the EPL distribution 
with its sub-models for a given data set. In fact to test w = 1, the LR statistic A = 
2{In L(a, B,o) —InL(a, B, 1)], where â, 8, and ô are the unrestricted estimates, and & and 
B are the restricted estimates. The LR test rejects the null hypothesis Ho if À > Xin where 
X; denotes the upper 100% point of the x? distribution with 1 degree of freedom. 
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7. APPLICATIONS 


In this section, the EPL distribution is applied to real data in order to illustrate the useful- 
ness and applicability of the model. We fit the density functions of the exponentiated power 
Lindley (EPL),power Lindley (PL), exponentiated Lindley (EL) and Lindley (L) distribu- 
tions. These models are also compared to the exponentiated Weibull (EW) and Weibull (W) 
distributions. Estimates of the parameters of the distributions, standard errors (in parenthe- 
ses), Akaike Information Criterion (AIC), Consistent Akaike Information Criterion (AICC), 
Bayesian Information Criterion (BIC) are given in Table 1 for the first data set, and Table 
2 for the second data set. 

The first data set is a subset of the breast feeding study from the National Longitudinal 
Survey of Youth, the complete data set is available in Klein, J.P., Moeschberger, M.L., 
Survival Analysis: Techniques for Censoring and Truncated Data, 2nd ed., Springer-Verlag 
New York, Inc., New York (2003). The data set considered here consists of the times to 
weaning 927 children of white-race mothers who choose to breast feed their children. The 
duration of the breast feeding was measured in weeks. 

The second data set is the Cancer Patients data. The data represents an uncensored data 
set corresponding to remission times (in months) of a random sample of 128 bladder cancer 
patients reported in Lee and Wang, (2003). It consists of the observations listed below: 0.08, 
2.09, 3.48, 4.87, 6.94, 8.66, 13.11, 23.63, 0.20, 2.23, 3.52, 4.98, 6.97, 9.02, 13.29, 0.40, 2.26, 
3.57, 5.06, 7.09, 9.22, 13.80, 25.74, 0.50, 2.46, 3.64, 5.09, 7.26, 9.47, 14.24, 25.82, 0.51, 2.54, 
3.70, 5.17, 7.28, 9.74, 14.76, 26.31, 0.81, 2.62, 3.82, 5.32, 7.32, 10.06, 14.77, 32.15, 2.64, 3.88, 
5.32, 7.39, 10.34, 14.83, 34.26, 0.90, 2.69, 4.18, 5.34, 7.59, 10.66, 15.96, 36.66, 1.05, 2.69, 
4.23, 5.41, 7.62, 10.75, 16.62, 43.01, 1.19, 2.75, 4.26, 5.41, 7.63, 17.12, 46.12, 1.26, 2.83, 4.33, 
5.49, 7.66, 11.25, 17.14, 79.05, 1.35, 2.87, 5.62, 7.87, 11.64, 17.36, 1.40, 3.02, 4.34, 5.71, 7.93, 
11.79, 18.10, 1.46, 4.40, 5.85, 8.26, 11.98, 19.13, 1.76, 3.25, 4.50, 6.25, 8.37, 12.02, 2.02, 3.31, 
4.51, 6.54, 8.53, 12:03, 20.28; 2.02, 3.36, 6.76, 12.07, 21.73, 2.07, 3.36, 6.93, 8.65, 12.63, 22.69. 

Estimates of the parameters of EPL distribution (standard error in parentheses), Akaike 
Information Criterion (AIC), Consistent Akaike Information Criterion (AICC) and Bayesian 
Information Criterion (BIC) are given in Table 8.1 for the first data set and in Table 8.2 for 
the second data set. 

Plots of the fitted densities and the histogram of the data are given in Figure 8.1 for the 
breastfeeding data, and Figure 8.2 for the cancer patients data. 

For the breasfeeding data, the LR statistic of the hypothesis Hp: PL(a,(,1) against 
Hy: EPL(a,8,w), is A = 7021.2 — 6981.3 = 39.9. The p-value is 2.67 x 1071° < 0.001. 
Therefore, we reject Ho in favor of Ha. Thus the EPL distribution performs better than the 
PL distribution. There is a significant difference between EPL and EL distributions with 
à = 33.9 and p-value=5.80 x 107° < 0.001. Thus, reject Hp: EL vs H, :EPL in favor of 
H,. A test of Ho: L vs H, : EL shows that À = 199.5 and p-value=2.68 x 1074. Thus, we 
reject Ho in favor of H,. To test the hypothesis Hp: EPL(a,3,w) against Ha: EW (a, p, w), 
we have A = 6981.3 — 6978.4 = 2.9. The p-value is 0.089 > 0.05. Thus we fail to reject Ho, 
and conclude that there is no significant difference between the EPL and EW distributions. 
There is a significant difference between EPL and W with the LR statistic A = 30.9 and 
p-value=2.7 x 1078. Thus the EPL distribution is preferred to the Weibull distribution. The 
values of the statistics AIC, AICC and BIC show that the EPL distribution is a “better” fit 
for the breast-feeding data. 
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For the Cancer Patients data set, the LR statistic for the hypothesis Hp: PL(a, 6,1) 
against Ha: EPL(a,6,w), is A = 5.8. The p-value is 0.016 < 0.05. Therefore, there is a sig- 
nificant difference between PL and EPL distributions. There is also a significant difference 
between EPL and EL distributions where À = 11.7. with a p-value of 6.25 x 1074 < 0.001. 
There is a significant difference between EL and L distributions with A = 6.5 and p- 
value=0.011 < 0.05. A test of Hp: EPL(a,8,w) vs Ha : EW(a,6,w) shows that \ = 0.5, 
and p-value=0.480 > 0.05, that there is no significant difference between the two distribu- 
tions. Based on the values of —2log L the EPL distribution fits the cancer patients data 
the best. However, the values of the statistics AIC, AICC and BIC are smaller for the EPL 
distribution and show that the EPL distribution is a “better” fit for the Cancer Patients 
data. 


8. CONCLUDING REMARKS 


We have presented and developed the mathematical properties of a new class of distribu- 
tions called the Exponentiated Power Lindley (EPL) distribution including the hazard and 
reverse hazard functions, moments, conditional moments, entropies, mean deviations, Lorenz 
and Bonferroni curves, Fisher information and maximum likelihood estimates. Applications 
of the proposed model to real data in order to demonstrate the usefulness and applicability 
of the class of distributions are also presented. 


TABLE 8.1. Estimates of Models for Breastfeeding Data 
Estimates Statistics 
Model a B Ww —2logL AIC AICC BIC 
EPL(a,6,w) 0.4186 1.0830 4.4536 6981.3 6987.3 6987.3 7001.8 
(0.0430) (0.2180) (1.3189) 

















PL(a, 8,1) 0.7349 0.2565 1 7021.2 7025.2 7025.2 7034.9 
(0.0158) (0.0128) 
EL, 8, w) 1 0.0837 0.5579 7073.4 7077.4 70774 7087.0 
(0.0033) (0.0242) 
L(1, B, 1) 1 0.1172 i 7272.9 7274.9 7274.9 7279.7 
(0.0027) 





EW(a,8,w) 0.5099 0.5625 3.7864 69784 69844 69844 6998.9 
(0.062) (0.1687) (1.1000) 

W(a,8,1) 0.9610 0.07017 1 7012.2 7016.2 7016.2 7025.8 
(0.0241) (0.0059) 








Density 
0.02 0.03 004 005 0.06 0.07 


0.00 0.01 
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Fitted Densities for Breastfeeding Data 








FIGURE 8.1. Plot of the fitted densities for the Breastfeeding Data 


TABLE 8.2. Estimates of Models for Cancer Patients Data 























Estimates Statistics 
Model a B w —2logL AIC AICC BIC 
EPL(a, 6,w) 0.5663 0.8191 2.7684 820.9 826.9 827.1 835.4 
(0.1017) (0.3116) (1.2903) 
PL(a, 6,1) 0.8302 0.2943 1 826.7 830.7 830.8 836.4 
(0.0472) (0.0370) 
EL(1, 6,w) 1 0.1649 0.7336 832.6 836.6 836.7 842.3 
(0.0166) (0.0912) 
L(1, 6, 1) 1 0.1960 1 839.1 841.1 841.1 843.9 
(0.0123) 
EW(a,8,w) 0.6544 0.4537 2.7960 821.4 827.4 827.6 835.9 
(0.1347) (0.2399) (1.2635) 
W(a, 6, 1) 1.0478 0.0939 1 828.2 832.2 832.3 837.9 
(0.0676) (0.0191) 
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Fitted Densities for Cancer Patients Data 
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FIGURE 8.2. Plot of the fitted densities for Cancer Patients Data 
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